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SECTION  1 


INTRODUCTION 

When  a  cable  lying  on  the  surface  of  the  earth  Is  exposed  to 
an  electromagnetic  pulse  from  a  nuclear  burst,  currents  are  Induced  on 
the  cable.  These  currents  can  be  on  the  order  of  a  few  hundred  to  a  few 
thousand  amperes,  depending  upon  such  things  as  soil  conductivity,  cable 
dimensions,  terminations,  and  angle  of  Incidence.  Surface  cables  are  used 
to  connect  various  types  of  tactical  equipment,  and  cable  responses  of 
this  magnitude  may  possibly  upset  or  damage  electronic  circuits  in  these 
equipments. 

Linear  analysis  £l-4]  of  surface  cable  response  to  both  high 
altitude  and  tactical  EMP  environments  predict  voltage  levels  that  may  be 
sufficient  to  break  down  cable  dielectrics.  In  addition,  the  predicted 
potential  gradient  across  the  soil  surface  Is  large  enough  to  cause  soil 
breakdown.  The  significance  of  the  effect  nonlinear! ties  have  on  the 
current  coupled  to  surface  cables  must  be  determined  and.  If  found  to  be 
significant,  must  be  Incorporated  Into  the  coupling  codes  used  for  surface 
system  EMP  studies. 

Nonlinear  cable  analysis  tools  have  recently  been  developed  for 
burled  cables  (5-7}.  One  such  tool  Is  the  computer  code  BLINE  which  was 
written  and  used  to  predict  these  effects.  The  BLINE  results  for  burled 
cables  have  shown  that  the  breakdown  of  soil  and  cable  dielectric  can  have 
an  effect  on  the  Induced  cable  currents.  The  general  Intent  of  this  study, 
therefore.  Is  to  do  a  similar  Investigation  Into  the  effects  of  nonlinearities 
on  surface  cable  response. 


Specifically,  the  study  has  the  following  three  objectives: 

(1)  Identify  and  document  all  nonlinear  parameters  associated 
with  the  coupling  of  EMP  energy  to  Army  system- type 
cables  when  they  are  exposed  to  "close-in"  tactical 
nuclear  environments. 

(2)  Determine  and  document  the  Importance  of  these  nonlinear 
parameters  in  terms  of  their  significance  to  the  magnitude 
of  energy  delivered  to  system  electronics  for  "close-in" 
EMP  coupling  situations. 

(3)  Consider  and  document  simple  experiments  that  may  be  used 
to  illustrate  the  nonlinear  effects  of  voltage  on  cable 
and  soil  parameters  used  in  calculating  the  coupling  of 
EMP  fields  to  surface  cables  exposed  to  "close-in" 
tactical  EMP  environments. 

In  view  of  the  above  objectives,  four  different  tasks  were 

addressed: 

(1)  Modification  of  BLINE  to  compute  surface  cable  response. 

(2)  Computation  of  linear  and  nonlinear  responses  for 
specific  cables. 

(3)  Determination  of  the  need  for  treatment  of  nonlinear 
air  conductivity. 

(4)  Determination  of  the  feasibility  of  using  high  level  EMP 
simulators,  such  as  AESOP  or  TEMPS  to  demonstrate  the 
expected  nonlinear  response  of  standard  cables. 

In  the  remainder  of  the  report,  modifications  of  BLINE  will 
first  be  addressed.  Then  the  computational  results  for  the  linear  and 
nonlinear  cases  will  be  discussed.  Then  the  feasibility  of  certain 
experiments  will  be  demonstrated.  Finally,  conclusions  will  be  made 
along  with  recommendations  for  further  investigations. 
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SECTION  2 


MODIFICATIONS  OF  BLINE  TO  COMPUTE  SURFACE  CABLE  RESPONSE 
2.1  Background 

BLINE  is  a  computer  code  developed  by  EMA  personnel  which 
solves  the  nonlinear  response  of  cables  buried  in  the  earth.  A  complete 
description  of  BLINE  can  be  found  in  Reference  5,  but  a  summary  will  be 
given  here  to  provide  enough  background  in  order  to  help  the  reader  under¬ 
stand  the  modifications  necessary  to  make  it  a  surface  cable  code. 

There  are  several  analysis  tools  used  to  predict  the  response 
of  buried  cables  to  EMP.  All  of  the  models  used  reduce  the  problem  to 
an  equivalent  transmission  line  (Figure  2.1).  In  the  equivalent  trans¬ 
mission  line  model  the  current  I  flowing  on  the  metal  conductor  is 
related  to  the  voltage  V  between  the  conductor  and  a  distant  reference 
point  by  the  telegraphers'  equations: 

-  -Z(u>)I(z,u»)  +  E^nc  (z, u>)  (2.1) 

and 

»  -Y(u>)  V(z,w)  +  Is(z,u)).  (2.2) 

Here  w  is  the  fourler  transform  variable,  Z(w)  and  Y(u>)  are  the  incre¬ 
mental  impedance  and  admittance  of  the  equivalent  line,  and  E1nc  and  I$ 
are  the  incident  electric  fields  and  impressed  current  density  (as  from 
a  lightning  stroke).  Because  these  parameters  account  for  the  presence 
of  both  the  cable  dielectric  and  the  soil,  any  simple  equivalent  (Figure 
2.1b)  must  have  elements  that  are  complex  functions  of  frequency.  For  the 
more  accurate  frequency  domain  treatment,  the  incremental  parameters  are 
determined  from  the  characteristic  Impedance  and  propagation  constant  of 

the  TM„  mode  as 
o 


a.  Original  Buried  Cable  Problem 


b.  Two-Equation  Transmission  Line  Equivalent 


Three-Equation  Transmission  Line  Equivalent 
Figure  2.1  Modeling  Buried  Conductor  Problem 


where  the  propagation  constant  y  and  characteristic  impedance  Z£  is 
solved  by  finding  zeros  of  a  complicated  transcendental  equation  [3]. 

For  more  approximate  frequency  domain  solutions  and  time  domain  approaches, 
the  use  of  three  telegraphers'  equations  is  more  appropriate  [8]: 


^  L—S’  -  “Z(w)I  (z»tu)  +  E1(IC  (z,«) 


|i  =  -Y j(ta))  Vj(z»w)  +  I.(z.u) 


=  -Ys  (to)  V$(z,u)  +  Is(z,a)). 


(2.4) 


(2.5) 


(2.6) 


In  this  formulation  the  parameters  Yd  and  Y$  are  more  easily  determined 
from  known  solutions  to  simpler  problems: 


Yd  (u)  d«f  Gd  +  jo)Cd 


+  ju 


'* ' '  ■  d  J  d 

Z(u)  dSf  Rc  +  juL  *  Rc  +  ju  in(^  z/\)  +  ^  *//a^ 


,  %  def 

Ys(u)  *  Ss  ♦  j«Cs 


+  ju 


r  a 
£n:  a  j 


(2.7) 


(2.8) 


(2.9) 


R  Is  the  DC  resistance  of  the  cable  conductor, 
c 

a  is  the  radius  of  the  cable  conductor, 

1 

a  Is  the  outer  radius  of  the  dielectric  sheath, 

2 

«s  Is  a  third  characteristic  dimension  that  depends  upon  the 
frequency  of  excitation. 


*VS  s 


.  <.  T.'  .  .  . 


Following  Vance  [93  one  would  use 

.  .  V?5  . 


r75T 


=  .794  6 


(2.10) 


The  time  domain  equivalent  of  the  telegraphers1  equations 
for  the  three  equation  models  are: 


3(Vd+V  )  3I  , 

— y~  *-lIt-rci+e 


(2,  t) 


aT  SVH 


(2.11) 


(2.12) 


dz  ■  “CS  at  "  Gsvs  + 


(2.13) 


Here,  the  values  of  Rc,  L,  Cd,  C$,  and  G$  are  taken  from  equations 
2.7,  2.8,  and  2.9  where  the  outer  radius  parameter  a^  is  not  time 
dependent.  It  is  noted  that  selecting  any  constant  value  overestimates 
the  inductance  (underestimates  the  induced  current)  at  early  time  (t<td) 
and  underestimates  the  Inductance  (overestimates  the  induced  current)  at 
late  times  (t<td),  where  td  is  determined  from  a  using: 

S  ■  rk  *  rrer^vir'  t2-14) 


m 

i  A+A 


2.2  Modifications  to  BLINE  to  Handle  Surface  Cables 

Because  BLINE  was  written  for  cables  completely  buried  in  the 
earth,  it  must  be  modified  to  treat  cables  lying  on  the  earth's  surface. 
Both  linear  and  nonlinear  aspects  must  be  taken  into  account. 

Modifications  for  the  linear  aspect,  which  for  our  purposes 
will  Include  the  effects  of  time  varying  air  conductivity,  involve  modi¬ 
fication  of  the  per  unit  length  Impedance  and  admittance  parameters.  For 
a  cable  burled  in  soil,  the  parameters  may  be  written 

zs  a  Rc  +  jwLs 

(2.15) 

Ys  *  «s  +  **f 
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In  the  above,  the  subscript  s  refers  to  the  soil  medium,  and  6S  is  the 
skin  depth  in  the  soil  given  by 


6  (a,)  =  V — £ — 
s'  ‘  wu  o 
o  s 


(2.17) 


In  order  to  convert  to  the  time  domain,  the  approximation  is  made  that 
t  *  2/f  [5],  so  that 


«s(t)  s 


V  SmT^- 


(2.18) 


Here,  t  is  measured  from  the  time  that  an  excitation  arrives  at  the 
location  of  Interest. 


In  an  Infinite  time  varying  conducting  air  medium  having  con¬ 
ductivity  aA(t),  formulas  analagous  to  those  for  soil  above  can  be  obtained 
merely  by  replacing  the  subscript  s  by  A. 

The  final  result  for  a  cable  on  the  surface  Is  then  the  combina¬ 
tion  of  the  above  results: 


surface 


zA(t)+  zs(t) 


(2.19) 


^surface  "  7  <V‘>+  V*)) 


These  formulas  then  define  the  equivalent  surface  cable  per  unit  length 
Impedance  and  admittance  including  the  effects  of  time  varying  air  con¬ 
ductivity.  They  can  also  be  applied  to  the  case  o.»o,  for  which  the  limit 


surface 


2Z$(t) 


(2.20: 


surface  "7s 


Y*(t) 


It  Is  noted  that  this  results  In  the  propagation  constant  of  the  surface 
cable  (In  the  presence  of  nonconducting  air)  being  equal  to  that  of  a 
burled  cable,  and  the  characteristic  Impedance  of  the  surface  cable  being 
twice  that  of  a  burled  cable. 

In  addition  to  modifying  the  soil  parameters  for  the  linear 
case,  the  dielectric  admittance  must  also  be  modified.  For  the  burled 
cable,  the  dielectric  capacitance  Is  simply  that  of  a  coaxial  cable 
formed  by  the  soil  and  conductor  boundaries.  The  same  Is  true  if  the 
cable  Is  at  the  Interface  between  two  conducting  media  such  as  conducting 
air  and  soil.  However,  as  the  air  conductivity  approaches  zero,  the  radial 
electric  field  lines  on  the  conductor  all  terminate  on  the  soil  and  are 
not  azimuthal ly  uniform  about  the  conductor. 


77777/7///' 


Figure  2.2  Geometry  Used  to  Calculate  Figure  2.3  Geometry  Used  to 


Soil  Per  unit  Length  Im¬ 
pedances  and  Admittances, 
and  Cable  Dielectric  Admit¬ 
tance  for  Conducting  Air. 


Calculate  Cable 
Dielectric  Admittance 
for  HAB  Case. 


Thus,  for  the  time  varying  conducting  air  medium  above  the  soil, 
the  geometry  of  Figure  2.2  can  be  used  to  compute  all  cable  parameters. 


--  -  -r * 


*  .  -  _  -  .  .  /.  ,'.  .\  .v  .V/ 


However,  for  no  air  conductivity  (the  HAB  case),  the  situation 
Is  somewhat  different.  A  surface  cable  would  not  be  deployed  In  the 
manner  of  Figure  2.2,  although  that  Is  the  correct  model  for  source 
region  computations.  (In  the  presence  of  conducting  air,  the  radial 
field  lines  should  be  azlmuthally  uniform,  even  if  deployed  like  Figure 
2.3.)  A  surface  cable  would  be  more  nearly  deployed  in  the  manner  of 
Figure  2.3.  It  Is  clear  that  In  this  case,  the  dielectric  capacitance 
Is  determined  by  the  geometry  of  a  wire  over  a  ground  plane.  This  also 
means  that  the  radial  electric  field  Is  not  simply  that  given  by  coaxial 
geometry  as  In  the  BLINE  case,  but  Is  given  by  the  wire  over  a  ground 
plane  geometry.  Figure  2.3.  The  electric  field  Is  therefore  a  maximum 
on  the  soil  side  of  the  conductor,  and  minimum  on  the  top. 

The  nonlinear  aspect  of  BLINE  to  be  modified  for  surface  cables 
Involves  the  treatment  of  soil  breakdown,  whereas  the  dielectric  break¬ 
down  Is  unmodified,  that  Is,  dielectric  breakdown  Is  treated  as  a  hard 
short. 

For  a  cable  burled  In  the  soil,  the  Ionized  soil  region  Is 
roughly  spherical  (except  for  a  bicone  coaxial  with  the  cable)  and  the 
radius  Is  determined- by  the  soil  breakdown  electric  field,  which  Is  l-2MV/m. 
For  a  cable  on  the  surface  of  the  earth,  an  arc  to  the  earth  will  form  a 
region  roughly  the  shape  of  a  thin  disc,  because  the  surface  breakdown  of 
soil  (250-500  kV/m)  Is  much  less  than  that  of  the  bulk  soil  medium.  There¬ 
fore,  the  Ionized  region  tends  to  extend  along  the  earth's  surface  Instead 
of  Into  It  [10], 

The  Ionized  region  Is  treated  as  a  thin  disc  of  radius 

arc’ 

where  r  Increases  until  the  electric  field  at  the  edge  of  the  disc 

flrC 

Is  equal  to  the  surface  breakdown  field  Egg.  For  a  current  I  flowing 
Into  the  disc,  r.„.  Is  determined  by 

orC 

Err  "  XJ  p  (2.21) 

BR  %  rarc 
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from  which  one  obtains 


arc 


\far 


s  *BR 


(2.22) 


given  by 


The  disc  then  has  a  conductance  and  capacitance  C__. 

ulC  Q(C 


arc 


4a  r 
s  arc 


arc 


4s  -r 
s  arc 


>  (2.23) 


The  relationship  of  the  conducting  disc  to  the  TEH  transmission 
line  surface  cable  code  is  Illustrated  in  Figure  2.4.  In  order  to  ensure 
that  r.^  is  never  greater  than  one-half  the  spatial  increment  Az,  r_„_ 

QiC  diC 

Is  set  equal  to  the  smaller  of  Az/2  or  r._. 

QiC 

It  is  noted  that  the  soil  may  break  down  before  the  cable 
dielectric  does.  When  the  dielectric  punctures,  current  flows  through 
the  puncture  into  the  soil  forming  a  disc  as  previously  described.  How¬ 
ever,  when  the  soil  breaks  down  and  the  dielectric  does  not,  the  concept 
of  the  disc  does  not  apply  because  the  current  flowing  into  the  soil  is 
distributed  along  the  cable  length  and  does  not  flow  from  a  discrete 
point.  In  this  case,  the  soil  nonlinearity  Is  treated  differently  as 
discussed  in  the  following. 


The  cable  geometry  of  Figure  2.2  Is  assured.  The  ionized  soil 
region  win  then  extend  into  the  soil  a  radius  of  a^  until  the  electric 
field  Is  equal  to  the  bulk  breakdown  field  of  the  soil.  This  breakdown 
field  Is  then  given  by 


Vso11 

a  An  S 
2  « 


(2.24) 


SECTION  3 


RESULTS  OF  SURFACE  CABLE  COMPUTATIONS 


3.1  Background 

One  of  the  goals  of  the  study  was  to  perform  surface  cable 
computations  for  several  cases,  with  and  without  Including  the  nonlinear 
effects.  Table  3.1  summarizes  the  cables  studied,  the  soil  conductivities, 
and  the  environments  used  for  the  parameter  study.  Both  linear  and  non* 
linear  predictions  were  made  In  order  to  ascertain  the  effects  of  Including 
non linearities. 


Table  3.1  Surface  Cables  for  Parameter  Study 


Length 

Envl ronment 

Soil  Conductivity  (mho/m) 

100  m 

HAB  Overhead 

.005 

HAB  Overhead 

.01 

Tactical  EMP 

.005 

Tactical  EMP 

.01 

600  m 

HAB  Overhead 

.005 

HAB  Overhead 

.01 

Tactical  EMP 

.005 

Tactical  EMP 

.01 

2800  m 

HAB  Overhead 

.005 

HAB  Overhead 

.01 

Direct  Hit  at  End  of  Cable 

.005 

Direct  Hit  at  End  of  Cable 

.01 

The  surface  cables  studied  are  formed  by  the  shield  and  outer 
dielectric  jacket  of  RG  58  and  RG  8  cables,  whose  dimensions  are  Indicated 
In  Table  3.2. 


r.|r 


111 


/ 

6 

2 

4.95 

9 

10.29 
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a.  wire  over  a  ground  plane 


b.  wire  in  conducting  medium 


Figure  3.1  Geometry  for  Definitions  of  Radial  Electric 
Fields. 

Table  3.3  Constants  Used  to  Multiply  the  Dielectric 
Voltage  to  Obtain  the  Radial  Electric  Field 
at  the  Points  Defined  In  Figure  3.1 


Cable 

RG  58 
RG  8 


The  radial  electric  fields  are  computed  by  a  constant  which 
multiplies  the  dielectric  voltages.  Figure  3.1  and  Table  3.3  define  the 
constants  at  the  locations  Indicated.  When  this  constant  is  multiplied 
by  the  dielectric  voltage,  the  result  Is  the  electric  field  at  that 
location  In  volts /meter. 

The  cable  jacket  material  Is  polyvinyl-chloride  (PVC).  In 
order  to  do  the  nonlinear  computations,  it  was  necessary  to  determine 
what  the  dielectric  strength  of  PVC  is.  Several  values  were  found  in 
the  literature  and  are  summarized  In  Table  3.4.  It  Is  clear  that  values 
vary  widely  and  depend  upon  the  types  of  material,  thickness  tested,  and 
test  waveform.  Values  for  polyethylene  (PE)  are  Included  for  reference 


purposes.  It  is  troublesome  that  the  relative  dielectric  strength  of 
PE  and  PVC  varies  considerably  from  source  to  source.  The  problem  is 
compounded  by  the  fact  that  for  the  purposes  at  hand,  it  is  desireable 
to  have  a  value  which  applies  to  a  voltage  which  rises  in  a  few  hundred 
nanoseconds  and  is  gone  nominally  within  one  or  several  microseconds,  arid 
probably  none  of  the  data  in  Table  3.4  applies  to  this  time  regime. 

Table  3.4  Summary  of  PVC  and  PE  Dielectric  Strength 

Data  from  Various  Sources.  Units  are  volts /mil. 


Source 

PVC 

PE 

Reference  12 

375 

420-550 

Reference  13 

400 

Reference  14 

600-1000  (flexible) 
700-1300  (rigid) 

460 

Reference  15 

800-1000 

460 

Reference  16 

800  (AC) 

1200  (AC) 

2000  (DC) 

3000  (DC) 

Extensive  data  has  been  taken  on  PE  breakdown  and  is  summarized 
in  Reference  17.  This  data  shows  that  for  lightning  type  pulses  (1*5  x 
40  ysec),  a  conservative  value  to  use  is  2000  V/mil.  It  should  be  noted 
that  some  data  exists  for  even  shorter  pulses  (1x5  usee)  which  yield  a 
5000  V/mil  value  [18J.  A  value  of  2000  V/mil  was  used  in  other  buried  cable 
studies  [5],  If  one  assumes  PVC  behaves  similarly  to  PE  and  uses  the  data 
from  Reference  16,  one  obtains  a  value  of  1330  volts/mil  for  PVC.  It  is 
admitted  that  this  is  a  guess  based  on  conflicting  data,  and  if  one  is 
going  to  perform  some  experiments  involving  cable  breakdown,  one  should 
first  obtain  the  dielectric  strength  data  for  the  time  regime  of  interest. 

It  is  noted  that  cable  manufacturers  were  called,  but  they  have  no  data  on 
the  jacket  dielectric  strength. 


3.2 


Computations  for  the  High  Altitude  Burst  (HAB) 


The  threat  for  the  HAB  predictions  was  taken  to  be  a  plane  wave 
normally  incident  upon  the  earth  and  polarized  parallel  to  the  cable. 
Previous  work  which  investigated  angle  of  incidence  effects  showed  that 
end-on  incidence  results  in  only  one-half  the  response  of  normal  incidence 
[3],  and  responses  for  other  angles  of  incidence  yields  values  within  one 
percent  of  the  normal  incident  case,  so  it  is  felt  that  normal  incidence 
is  the  appropriate  case  to  study.  The  temporal  behavior  of  the  plane 
wave  is  the  familiar  double  exponential 


^inc 


=  52  (e 


-4x10  H  .  ,-4.76x10*1,  kV/m 


(3.1) 


Standard  plane  wave  reflection  coefficients  were  used  to  obtain  the  total 
field  at  the  earth's  surface,  which  is  shown  in  Figure  3.2  for  the  two 
values  of  conductivity  used  in  the  study.  The  cables  studied  were  shorted 
on  one  end  and  open  on  the  other. 


Figure  3.2  Total  HAB  Electric  Field  at  Earth's  Surface 
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HA8  results  are  summarized  in  Figures  3.3  through  3.8.  Some 
general  remarks  can  be  made  about  the  HAB  results.  First,  it  was  found 
that  the  results  for  lengths  100  m  or  greater  were  independent  of  length, 
within  1%.  That  is,  the  spectral  content  of  the  HAB  input  is  of  suffi¬ 
ciently  high  frequency  so  that  signals  propagating  from  one  end  toward 
the  other  attenuate  rapidly  and  are  not  discernible  at  the  far  end.  This 
is  evident  in  Figure  3.3  which  shows  the  short  circuit  current.  The 
response  from  the  other  end  should  be  arriving  at  about  700-800  ns  assuming 
a  nominal  propagation  velocity  of  one-half  that  of  light,  and  it  is  clear 
that  no  such  signal  is  there. 

The  second  observation  is  that  the  voltage  induced  on  the  R6  8 
dielectric  is  insufficient  to  cause  breakdown  (52  kV  is  required),  and 
RG  58  breaks  down  only  for  the  case  a  =  .005  mho/m. 

From  Figure  3.3  it  is  evident  that  the  nonlinearities  (both 
soil  and  dielectric)  do  not  make  a  significant  difference  in  the  short 
circuit  current  response.  For  this  case,  only  the  three  cells  closest 
to  the  open  end  broke  down,  and  this  effect  is  not  felt  at  the  shorted 
end  nearly  100  m  away.  The  three  cells  closest  to  the  open  end  broke 
down  at  267,  172,  and  270  ns,  ordered  in  the  direction  from  the  shorted 
end  to  the  open  end,  respectively. 

It  is  also  evident  that  the  soil  nonlinearities  (which  occur 
even  if  the  dielectric  does  not  break  down)  are  unimportant  in  the 
current  response  as  evidenced  in  Figures  3.3  and  3.6.  It  also  makes 
little  difference  in  the  soil  and  dielectric  voltages  as  evidenced  in 
the  RG  8  responses  of  Figures  3.7  and  3.8,  which  show  essentially  the 
same  response  for  linear  and  soil  nonlinear  cases. 

For  the  case  in  which  the  RG  58  dielectric  breaks  down,  the 
voltages  across  the  soil  and  dielectric  are  quite  different  from  the 
linear  case  as  one  would  expect,  and  this  difference  is  evidenced  in 
Figures  3.4  and  3.5.  It  is  noted  the  end  of  the  cable  was  artificially 
forced  to  be  open,  which  is  why  the  dielectric  voltage  of  Figure  3.5  does 
not  go  to  zero. 
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Finally,  it  is  noted  that  the  voltages  induced  on  the  cable 
are  sufficient  to  cause  air  breakdown,  which  was  not  included  in  the 
model.  For  RG  58  from  Figure  3.5  and  Table  3.3,  one  gets  an  electric 
field  in  the  air  above  the  cable  more  than  7  MV/m,  which  is  sufficient 
to  cause  air  breakdown.  For  RG  8,  the  field  is  larger  than  3.6  MV/m, 
which  is  still  more  than  enough. 


o  =  .005,  Linear  and 
Nonlinear 


0.01,  Linear  and 
Nonlinear 


time,  ns 


Figure  3.3  HAB  Short  Circuit  Current  on  100  m  RG  58  Cable 


Figure  3.8  HAB  Dielectric  Voltage  on  Open  End  of  100  m  RG  8. 

Linear  and  Nonlinear. 

3.3  Computations  for  a  Tactical  Environment 

The  tactical  environment  which  excites  and  influences  the  cable 
response  Is  the  radial  electric  field  and  the  air  conductivity.  HDL  sup¬ 
plied  waveforms  for  these  variables  at  ranges  from  50  to  2800  meters  which 
were  used  in  the  computations.  Typical  waveforms  are  shown  in  Figures 
3.9  and  3.10  which  are  for  a  range  of  1000  m.  Cables  of  100,  600,  and 
2800  m  length  were  studied.  In  all  cases,  they  are  oriented  radially  to 
the  source,  and  the  end  closest  to  the  source  Is  a  perfect  short  and  the 
end  farthest  away  from  the  source  Is  open  circuited  (except  for  one  2800  m 
case,  discussed  later).  Responses  computed  are  the  short  circuit  current 
at  the  shorted  end,  and  the  dielectric  voltage  and  soil  voltage  at  the 
open  end. 
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Results  for  the  100  m  cable  are  shown  in  Figures  3.11  through 
3.16.  The  cable  is  located  between  the  ranges  of  900  and  1000  m.  The 
effect  of  the  nonlinearities  is  to  greatly  increase  the  shor*  circuit 
current,  as  shown  in  Figures  3.11  and  3.12.  The  late  time  i^ponse  is 
also  greatly  extended.  The  effect  of  nonlinearities  on  the  soil  voltage 
Is  to  reduce  It  and  significantly  alter  the  spectral  content  as  evidenced 
in  Figures  3.13  and  3.14.  Of  course,  as  one  would  expect,  the  effect  of 
the  nonlinearity  Is  to  greatly  decrease  the  dielectric  voltage. 

The  propagation  of  dielectric  breakdowns  for  the  100  m  cable 
is  Interesting.  For  example,  for  the  RG  58  at  a  •  .005,  the  dielectric 
breakdowns  begin  at  the  open  end  at  t  ■  1.63  usee  and  propagate  a  distance 
of  38.7  meters  toward  the  source,  at  which  point  breakdown  occurs  at  7.95 
usee.  It  appears  likely  that  If  the  computation  were  extended  in  time, 
the  punctures  would  propagate  even  further. 

Results  for  the  600  m  cables  are  shown  in  Figures  3.17  through 
3.22.  One  end  Is  at  a  range  of  400  m  and  the  other  is  at  1000  m.  Results 
are  similar  In  character  to  those  of  the  100  m  cables  except  that  the 
responses  are  much  larger.  The  propagation  of  dielectric  breakdowns  is 
more  complex  In  this  case  than  for  the  100  m  case.  For  RG  58  at  a  *  .005, 
the  first  puncture  occurs  at  92  meters  from  the  shorted  end  at  3.23  usee 
and  propagates  both  directions;  inward  to  74  meters  at  3.3  usee,  and  then 
outward  every  6.6  meters  (approximately).  In  addition,  breakdown  occurs 
at  the  far  end  at  3.37  usee  and  propagates  inward  a  distance  of  66  meters 
where  breakdown  occurs  at  12.6  usee.  The  end  result  is  a  complex  pattern 
of  punctures  which  is  non-uniformly  distributed  along  the  cable.  There 
are  no  punctures  at  distances  closer  than  74  meters  to  the  source  end  of 
the  cable. 


a  »  0.01 
Nonlinear 


a  ■  0.005 
Linear 


o  »  0.005 
Nonlinear 

o  ■  0.01 
Linear 


Figure  3.11  Tactical  Short  Circuit  Current  for  100  m 
RG  58  Cable. 


a  *  0.01  Nonlinear 


*  0.05 

^  Nonlinear 

a  ■  0.009\  a  »  0.01 
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Figure  3.12  Tactical  Short  Circuit  Current  for  100  m 
RG  8  Cable. 
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Figure  3.13  Tactical  Soil  Voltage  for 
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Figure  3.21  Tactical  Dielectric  Voltage  for  600  m  RG  58 


Figure  3.22  Tactical  Dielectric  Voltage  for  600  m  RG  8. 
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Results  for  the  2800  m  cable  are  given  in  Figures  3.23  through 
3.28.  Cable  currents  are  similar  in  character  and  in  amplitude  to  those 
for  the  600  m  cable.  The  effect  of  the  nonlinearities  is  to  greatly 
increase  the  late  time  response.  Also,  the  nonlinearities  greatly  reduce 
the  dielectric  voltages  and  significantly  alter  the  amplitude  and  wave¬ 
shape  of  the  soil  voltages. 

The  pattern  of  dielectric  punctures  is  quite  complex  for  this 
case  also.  For  the  RG  58  with  a  =  .005,  the  first  puncture  occurs  at 
280  m  at  3.18  usee.  The  puncture  closest  to  the  source  is  at  60  m  and 
occurs  at  3.3  usee.  The  far  end  breaks  down  at  13.8  usee,  but  this  is 
the  only  puncture  in  the  last  500  meters  of  the  cable.  The  remainder  of 
the  cable  has  a  complex  distribution  of  punctures. 

For  all  of  the  tactical  computations,  voltage  levels  are 
higher  than  those  obtained  for  the  HAB  case.  Thus  it  is  clear  that  the 
air  will  break  down  also,  an  effect  not  included  in  these  computations 
but  which  should  be  considered  in  the  future. 

Because  there  is  some  interest  in  knowing  what  the  short 
circuit  current  at  the  far  end  of  the  2800  m  cable  is,  it  was  also 
calculated  (with  both  ends  shorted)  and  is  shown  in  Figure  3.29.  It 
is  noted  that  the  linear  and  nonlinear  cases  agree.  This  is  because 
there  is  no  dielectric  breakdown  in  the  last  500  meters  of  the  cable, 
and  the  effects  of  breakdown  are  not  significant  by  the  time  they 
propagate  to  the  far  end. 


Figure  3.26  Tactical  Soil  Voltage  for  2800  m  RG  8  Cable 


SECTION  4 


EXPERIMENT  FEASIBILITY 


4.1  Background 

One  of  the  goals  of  this  study  was  to  determine  the  feasibility 
of  using  existing  simulators  to  investigate  nonlinear  cable  response. 
Computations  were  done  using  AESOP,  which  has  the  same  electromagnetic 
characteristics  as  TEMPS,  as  the  source  of  EMP  fields  to  illumine  the 
cables  into  their  nonlinear  response  regimes.  It  was  found  that  AESOP 
can  indeed  be  used  to  perform  these  experiments.  In  order  to  do  this 
study,  the  electromagnetic  field  model  of  AESOP  had  to  be  determined. 

Then  some  experimental  concepts  are  given,  and  then  some  computational 
results  are  shown. 

For  the  purposes  of  this  section,  three  cable  responses  were 
computed,  both  with  and  without  nonlinear  effects  included.  Figure  4.1 
illustrates  the  type  of  cable  configurations  that  were  studied.  A  small 
coaxial  cable  like  RG  58  would  make  a  reasonable  test  cable  for  the 
insulated  cable  and  1/2"  electrical  conduit  would  suffice  for  the  uninsu¬ 
lated  conduit. 

The  three  experiments  include  most  of  the  interesting  effects. 
The  U-shaped  insulated  cable  will  exhibit  a  voltage  maximum  at  the  base 
of  the  U.  If  breakdown  occurs,  it  will  begin  in  the  U-shaped  region. 
Several  punctures  of  the  dielectric  are  likely  to  occur.  The  principal 
question  is  at  what  minimum  field  intensity  will  a  puncture  of  the 
insulating  jacket  be  obtained.  Additionally,  one  would  like  to  know  the 
periodicity  of  breakdowns  along  the  cable  in  the  case  where  a  very  high 
field  intensity  pulse  is  applied  to  an  undamaged  cable. 


Cable  Studies 


The  uninsulated  conduit  will  be  a  test  of  our  understanding 
of  the  nonlinear  response  of  soil,  unconfused  by  the  mixing  of  dielectric 
breakdown  and  soil  breakdown. 

The  ground  stake  experiment  will  examine  the  ground  stake 
response  as  a  termination  Impedance.  This  tests  directly  the  model  of 
nonlinear  incremental  soil  parameters  without  including  distributed 
parameters  and  distributed  sources  in  the  analysis. 

4.2  The  AESOP  Field  Model 

The  model  used  to  predict  the  fields  is  the  same  as  that  of 
Reference  19.  The  electromagnetic  fields  produced  by  AESOP  are  much 
more  complicated  than  the  high  altitude  burst  plane  waves.  For  AESOP, 
the  amplitude,  angle  of  incidence,  and  the  wave  polarization  are  a 
function  of  position,  whereas  they  are  not  for  HAB. 

Figure  4.2  shows  an  oblique  view  of  AESOP  and  the  coordinate 
system.  The  antenna  Is  on  the  z  axis  at  y  *  x  *  0.  It  Is  20  meters  above 
the  earth.  If  one  draws  a  vector  from  the  origin  at  (0,  0,  0)  to  another 
point  (x,  y,  z)  where  the  fields  are  to  be  defined,  it  is  useful  to  define 
the  angles  B  and  $  according  to  the  following: 


The  radiated  fields  are  found  from  the  formulas  of  Jordon  [20] 
for  fields  near  a  dipole  and  are  expanded  In  a  power  series  to  obtain 
the  fields  before  reflection  from  the  ends  take  place.  This  is  equivalent 
to  making  the  antenna  Infinitely  long,  or  else  having  a  long  antenna  with 
perfectly  matched  terminations  on  each  end. 


The  resultant  total  electric  and  magnetic  fields  are 

V 


Et  ■ 


•«ii  V  x2  +  y2 


(4.2) 


where  VQ  is  the  pulser  voltage  and  ^  Is  a  constant,  and  are  chosen  to 
normalize  the  incident  electric  field  to  50  kV/m  at  the  earth  at  the 
coordinate  y  =  50  m.  The  field  and  its  spectral  content  are  shown  in 
Figure  4.3  and  4.4,  respectively. 


given  by 


The  unit  vector  Up  perpendicular  to  the  plane  of  incidence  is 
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Then  the  incident  horizontally  polarized  and  vertically  polarized  inci¬ 
dent  fields  are  given  by 
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An  approximation  is  made  that  the  plane  wave  reflection  co¬ 
efficients  [21]  for  the  earth  can  be  used.  The  total  field  parallel  to 
the  surface  cables  is  then  computed  and  used  as  the  excitation  source 
which  drives  them.  The  earth  is  assumed  to  have  constitutive  parameters 
e  «  10^  and  a  *  .005  mho/m.  The  risetime  of  the  AESOP  field  can  be  in 
the  range  of  4-12  ns.,  and  a  nominal  value  of  8  ns  was  used  in  the  compu¬ 
tations. 


Figure  4.5  shows  the  total  electric  field  at  (-20,  3,  0).  It 
was  found  that  along  a  line  from  this  point  to  (-20,  3,  37),  the  total 
electric  field  did  not  noticeably  change  waveshape  and  the  peak  amplitude 
decreased  nearly  linearly  (within  3%)  to  a  minimum  of  45.2  kV/m.  This 
line  is  the  line  upon  which  all  of  the  cables  to  be  studied  were  placed. 
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Figure  4.5  AESOP  Total  Electric  Field  at  (-20,  3,  0). 


4.3  Results  for  the  Ground  Stake  Experiment 

The  ground  stake  experiment  consisted  of  37  meters  of  RG  58 
with  its  shield  terminated  to  a  1"  OD  ground  stake  of  length  2  meters 
and  the  other  end  open  circuited.  The  ground  stake  impedance  is  98  a, 
according  to  the  formulas  of  Sunde  [10],  Figure  4.6  shows  the  stake 
current  and  Figure  4.7  shows  the  dielectric  voltage  at  the  other  end  of 
the  cable  for  linear  and  nonlinear  cases. 

The  current  shows  a  small  nonlinear  effect,  but  it  would  be 
noticeable  on  a  measurement.  If  one  could  measure  the  dielectric  voltage 
a  marked  difference  could  be  observed. 


Figure  4.6  Stake  Current 
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Figure  4.7  Dielectric  Voltage  at  End  of  the  Insulated  Cable. 

The  computations  indicate  that  the  7  meters  of  cable  closest  to  the  open 
end  would  experience  multiple  punctures.  The  first  puncture  occurs  at 
the  open  end  at  84  ns  and  the  punctures  propagate  7  meters  and  the  last 
one  occurs  at  330  ns. 

Because  punctures  occur  only  at  the  end  away  from  the  stake, 
their  effect  on  the  stake  current  is  rather  small.  A  shorter  cable  would 
probably  cause  a  bigger  difference  in  stake  current,  but  even  that  in 
Figure  4.6  is  discernible.  From  Figure  4.7,  one  may  conclude  that  if  one 
reduced  the  simulator  output  by  a  factor  of  2,  the  nonlinearities  would 
just  be  at  the  threshold  of  occurrence. 


4.4 


Bare  Conduit  Experiment 


Figures  4.8  through  4.10  present  the  results  for  the  37  meter 
long  bare  conduit  computations.  These  results  are  consistent  with  pre¬ 
vious  results  in  that  only  small  effects  in  cable  response  are  caused 
by  soil  nonlinearities.  The  results  shown  are  computed  from  a  much 
reduced  value  of  bulk  soil  breakdown  voltage:  375  kV/m.  This  was  done 
because  if  one  used  2  MV/m,  the  results  nearly  overlayed  each  other  and 
it  was  desirous  to  ascertain  the  importance  of  the  soil  breakdown  field. 
It  appears  that  the  non-sensitivity  of  response  to  soil  breakdown  above 
can  be  readily  checked  by  this  experiment. 
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Figure  4.8  Conduit  Center  Current 
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4.5 


"U"  Cable  Experiment 

The  "U"  cable  is  37  meters  long  on  each  side  and  the  joining 
part  is  6  meters  long.  The  cable  is  open  on  the  ends.  Its  response  is 
predicted  to  show  the  most  dramatic  nonlinear  effects.  Results  are 
presented  in  Figures  4.11  through  4.16. 

The  current  response  (Figures  4.11  through  4.13)  should  be 
easily  measurable.  The  effect  of  the  nonlinearities  is  to  greatly  increase 
the  peak  value  and  late  time  currents  as  shown.  From  the  voltage  waveforms 
one  can  estimate  that  nonlinearities  would  begin  to  occur  at  about  50%  of 
simulator  output. 


Figure  4.11  Current  at  Corner  of  "U"  Cable. 


Figure  4.16  "U"  Cable  End  Dielectric  Voltage. 

For  the  MU“  cable,  the  punctures  are  in  the  vicinity  of  the 
corners  and  at  the  open  ends.  At  each  open  end,  a  puncture  occurs  at 
155  ns  and  punctures  propagate  away  from  the  end  a  distance  of  2.6  m, 

where  the  last  puncture  occurs  at  491  ns,  (time  measured  from  when  the 

» 

incident  wave  arrives  at  the  comer  of  the  "U").  In  the  cable  region 
between  the  two  comers  there  are  numerous  punctures  which  occur  between 
153  ns  and  177  ns.  It  is  interesting  to  note  that  no  punctures  propagate 
away  from  the  corners  towards  the  open  ends. 


SECTION  5 


CONCLUSIONS 

This  is  the  first  known  attempt  to  estimate  the  effects  of 
soil  and  dielectric  breakdown  upon  surface  cable  transient  response, 
and  the  results  need  to  be  subjected  to  careful  experimental  comparisons. 
However,  the  analytical  results  indicate  several  interesting  conclusions: 

1)  The  effect  of  iionlinearities  is  principally  to  greatly 
increase  the  late  time  currents  and  reduce  the  voltage  delivered  by  a 
cable  (for  the  tactical  case).  This  means  that  the  effect  of  nonline¬ 
arities  on  system  response  is  to  reduce  the  voltage  on  high  impedance 
systems  and  to  increase  the  current  on  low  impedance  systems.  Of  course, 
a  high  impedance  system  may  itself  break  down  and  become  a  low  impedance 
system. 

2)  HAB  results  were  not  greatly  affected  by  nonlinearities, 

at  least  for  the  cables  studied.  More  dramatic  effects  might  be  observed 
on  wires  with  thinner  Insulation,  such  as  field  wire. 

3)  The  most  Important  nonlinearity  is  cable  dielectric  break¬ 
down,  whereas  soil  breakdown  Is  of  minor,  and  sometimes  negligible, 
importance. 

4)  The  results  clearly  show  that  nonlinear  air  breakdown  will 
occur  in  the  vicinity  of  the  cable  for  both  HAB  and  tactical  excitations. 
This  effect  should  be  evaluated  as  a  next  step,  especially  before  performing 
experiments  In  TEMPS  or  AESOP. 

5)  One  of  the  most  Important  unknowns  is  the  dielectric  break¬ 
down  of  cable  insulations  In  the  appropriate  time  regime.  Data  from  the 
literature  and  handbooks  is  Incomplete  and  conflicting.  It  is  felt  that 
the  parties  concerned  should  take  their  own  data  in  order  to  have  confidence 
in  it. 
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6)  AESOP  or  TEMPS  can  be  used  to  validate  the  nonlinear 
analysis  tools  used  in  this  study.  Computations  show  that  the  important 
features  of  nonlinear  cable  response,  such  as  the  late  time  large  current 
and  the  relative  importance  of  the  soil/dielectric  breakdown,  can  be 
determined. 

7)  It  is  noted  that  results  were  presented  for  a  cable  "ideally" 
located  at  the  earth's  surface.  In  reality,  a  cable  will  perhaps  be 
touching  the  earth's  surface  at  only  randomly  selected  points,  and  will  be 
slightly  elevated  for  random  intervals.  Results  should  be  computed  for 
other  representative  geometries. 

It  should  be  pointed  out  that  there  are  several  assumptions  in 
the  model  which  may  require  further  research  or  verifications.  Among 
these  are: 

1.  The  value  for  the  soil  breakdown  fields. 

2.  The  temporal  behavior  of  soil  breakdown  (that  is,  it 
was  assumed  that  the  soil  broke  down  without  time 
delay,  once  the  field  exceeded  the  breakdown  value). 

3.  The  validity  of  the  "disc  model". 

4.  The  "disc"  radius  was  limited  to  be  one- half  of  a 
spatial  finite  difference  increment. 

5.  The  value  of  the  cable  dielectric  breakdown  fields 
and  their  temporal  dependence. 

6.  The  validity  of  the  model  used  for  soil  breakdown 
when  the  dielectric  did  not.  (Surface  flashover 
was  not  assumed.) 

7.  The  accuracy  of  replacing  f  by  2/t  in  the  skin  depth 
formulas  (Equation  2.18). 
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